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Abstract 

A generalized FRW model of space-time is studied, taking into considera- 
tion the anisotropic structure of fields which are depended on the position 
and the direction (velocity). The Raychaudhouri and Friedman- like equa- 
tions are investigated assuming the Finslerian character of space-time. A 
long range vector field of cosmological origin is considered in relation 
to a physical geometry where the Cartan connection has a fundamen- 
tal role. The Friedman equations are produced including extra anisotropic 
terms. The variation of anisotropy z t is expressed in terms of the Cartan 
torsion tensor of the Finslerian manifold. A physical generalization of the 
Hubble and other cosmological parameters arises as a direct consequence 
of the equations of motion. 
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1 Introduction 



During the last few years considerable studies concerning observable anisotropics 
of the universe have been investigated^ [7] [T8l [19] .These are connected to the 
very early state of the universe and related to the estimations of WMAP of CMB, 
the anisotropic pressure or the incorporation of a primordial vector field(e.g. a 
magnetic field) to the metrical spatial structure of the universe[91 EH 150] -In 
this case the form of scale factor can be influenced by the introductory field. A 
geometry which may connect the Riemannian metric structure of the space-time 
to physical vector fields, is the class of Finsler-Randers type spaces. In these 
spaces an electromagnetic field, a magnetic field or a gauge vector field may 
emerge out by a physical source of the universe and can be incorporated into 
the geometry causing an anisotropic structure [2131 |2"TI |2"21 |2"T1 |2"B1 I3T1 135] . 

Finsler geometry or the theory of Finsler spaces may be considered as a 
generalized Riemannian geometry of the first order within the sphere of metrical 
differential geometry [33] . A Finsler space is a metric space in which the metric 
function is defined by a norm T on a tangent bundle instead of defining an 
inner product structure on it. The norm will be a real function T(x,y) of a 
space-time point x and a tangent vector y £ T X M which places the role of an 
internal variable (Appendix A) . This y dependence characterizes essentially the 
Finslerian field and has been combined with the concept of anisotropy which 
causes the deviation from Remannian geometry [48] 149] .All kinds of generalized 
metric theories and unified field theories belong to the larger class of the so 
called anisotropic field theories 48, 49] . Therefore these geometrical anisotropies 
are caused by internal variables. Under these conditions a Finsler geometry can 
be considered as a physical geometry on which matter dynamics takes place 
while the Riemann geometry is the gravitational geometry [231 |2"T1 133] . 

The Cartan's torsion tensor [33J [35] characterizes all the geometrical concept 
of Finsler geometry and appears to all expressions of geometrical objects such 
as connection and curvature. In some cases it is useful from a physical point 
of view consider a vector field in the form y % {x)ji = 1,2,3,4 and the induced 
Finslerian metric tensor gives rise to the osculating Riemannian metric tensor 
r^ v (x) = fnu(x,y(x)). 

In the present paper we adopt such an approach in order to obtain some 
results concerning to a Friedman-like Robertson Walker cosmological model(F- 
LRW) and its consequences. We proceed by introducing a Randers-type La- 
grangian and the induced Finslerian metric modified appropriately for a weak 
primordial vector field u a 28\ (sec. 2). We construct the field equations using the 
osculating approach [33J [35J and derive a Friedman-like equation of motion with 
an extra anisotropic term (sec. 3, 4). We generalize the Hubble parameter H, the 
density parameter O and the deceleration parameter q for a weak anisotropic 
approach. The cosmological parameters depend directly on the anisotropy gen- 
erated by the vector field defined in the Lagrangian metric function (sec. 4). 
The solutions of the Fricdman-like equation for both matter and radiation dom- 
inated universe, the Raychaudhuri equation initially presented in [211 122] and 
the CMB temperature estimation are all affected by the presence of the rate of 
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anisotropy at the field equations (sec. 4, 5). The anisotropic solution of the scale 
factor coincides with the standard ones derived under the assumptions of homo- 
geneity and isotropy. A possible estimation of the rate of anisotropy parameter 
might be possible if we consider intermediate values for the cosmological con- 
stant for a generalized de-Sitter model defined by the Friedman-like equation 
of motion with a cosmological const ant |44). A more accurate estimation of the 
model's extra parameter can be made by calculating the CMB shift and the 
baryon oscillation acoustic peak using the Friedman equations and comparing 
them to the corresponding values coming from the WMAP data[2l 151 [SI IT21 [T5] . 
The data produced by these tests depicts a flat universe confirming the obser- 
vations. 



2 Anisotropy and a Randers type Lagrangian 
metric 

An alternative way of studying of physical phenomena is to incorporate the 
dynamics to the active geometrical background following Einstein's meaning 
of gravity. Our investigation is based on the introduction of a Lagrangian 
metric [221 [5S] considering the anisotropy of the universe [T] as an embodied char- 
acteristic of the geometry of space-time. A similar investigation has been applied 
to the case of electromagnetism[2ni OH] together with some recent progress in 
gravity, cosmology and fluid dynamics [2"T1 |2"51 |2"5] . We consider the geodesies of 
the 4-dimensional space-time to be produced by a Lagrangian identified to be 
the Randers-type metric function(the greek indices belong to {0, 1, 2, 3} and the 
latin ones to {1, 2, 3}) 

F(x,y)= a{x, y) + 4>{x)k a y a (2.1) 
v{x,y)= v / a KX (x)y K y x (2.2) 

where a K \(x) is the Robertson- Walker metric defined as 

a 2 

cl k \{x) — diag(l, — r, — a 2 r 2 , —a 2 r 2 sin 2 9). (2-3) 

1 — kr 

where k = 0, ±1 for a flat, closed and hyperbolic geometry respectively. The 
spatial coordinates are comoving and the time coordinate represents the proper 
time measured by the comoving observer. The vector y M = represents the 
tangent 4-velocity of a comoving observer along a preferred family of world- 
lines(fluid flow lines) in a locally anisotropic universe;the arclength parameter 
s stands for the proper time. We proceed by considering the natural Lorentzian 
units i.e. c = 1. If we fix the direction y = x then o~(x,x) — l.The vector field 

u a (x) — k a (f>(x) (2-4) 

stands for a weak primordial vector field \u a \ <C 1 incorporated to the geometry 
of space-time as an intrinsic characteristic. This field would most naturally be 
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expected to point in the same direction with the tangent vectors of the fluid 
flow lines [41]. As a result it will have only a time like component which can 
be expressed as a function of the proper time u a = (uq, 0,0,0). Important 
information about the anisotropy is encoded into the component Mn(t)|41j . We 
consider a linearized variation of anisotropy, therefore the approximation 

cj>{x) « ^(0) + 9 M ^(0)/ (2.5) 

is valid for small x. 



The choice of the connection A* (x) and the 
curvature. 



The metric of the Finsler space can be directly calculated from the metric func- 

1 d 2 F 2 
1 dy^dy" 



tion F. Since f^ v (x,y) = \ S'^ (x, y) we derive 



L(- 3 



ffiu = g^v + —( U v-Vv + u vVn) aVuVu + (3.1) 



where 



F 

g»v{x,y) = — (x,y)a^(x) (3.2) 

(7 

and 

0(x,y) = (j){x)k a y a = u a (x)y a (3.3) 

Under the weak field assumption we can approximate the Finslerian metric 
(|3.1|) as a perturbation of the FRW metric since in General Relativity a weak 
vector field in a space (e.g. primordial magnetic field) can be treated as first 
order perturbation of the Riemann metric tensor. The metric is considered 
to have signature (+,—,—,—) for any (x,y). The square of the length of an 
arbitrary contravariant vector is to be defined \X\ 2 — f^ u {x, y)X^X u . The 
connection components of the metric are given by (|A.7[) . 

In many cases we consider a convenient Finsler metric to approximate the 
gravitational theories [33J 13S3- This metric is connected to a Riemannian one, 
Tfj.u(x) referred as osculating Riemannian metric\33\ 

»>/(x) = f liv {x,y{x)) (3.4) 

with the following Christoffel components 

dv p 

r1A x ) = 7A M (z,y(z)) +C; p {x,y{x))-^ x (x) + 

+C^(x,y(x))^(x)-g^(x,y(x))C^ p (x,y(x))^(x) (3.5) 
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thus the equation of geodesies is given by 



^+^(^=0. (3.6) 

Under the assumption that the vector field y a satisfies the relation y? v = the 
Finslerian 8— covariant derivative and the Cartan's covariant derivative of an 
arbitrary vector field X a (x) are equal [33l [35] (see Appendix B) 

X%(x,y(x)) = Xfa(x,y(x)). (3.7) 

The Cartan's torsion tensor can be easily deduced from (|2.1[) and (|A.8p the full 
expression is [25] 



C^vx = \ J — S^ v x){a^ v ux) - -^S^ y x){y^yuux) - ■^S^x)(a^yx)\ (3.8) 

£ I (7 (7 (7 J 

where S^x denotes the sum over the cyclic permutation of the indices. Every 
single term of (|3.8p is proportional to the components of the field u a thus C m „a ~ 
under the condition \u a \ <C 1 and then we can drop all the torsion dependent 
terms in (13 . 5|) . Therefore the approximation for the Christoffel components 
becomes 

where A Xv {x) represent the osculating affine connection coefficients. The affine 
curvature tensor associated with the proper choice of the connection coefficients 
A\ v gives directly the curvature which is associated with the commutation re- 
lations of the ^-derivatives 

Lxpu = Ax u ^ — j4a M)W + A p Xu A K pil — A p Xp A K pv . (3.10) 

The Ricci tensor is given by 

= L a pav (3-11) 



and the scalar curvature 



L = r u L^. (3.12) 



The inverted metric is calculated in[5S]. The components of the Ricci tensor 
can be simplified due to the conditions 

iio-O (3.13) 

ill w 0. (3.14) 



The condition (|3.13[) is valid since <f>(x) can be written at the linear form (|2.5p 
together with (|3.14|) where we have considered u very small at the first stages 
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of a highly accelerated expanding universe [38]. We arrive then at the following 
nonzero components 




3 (a/a + 3/4d/awo) 

- (ad + 2d 2 + 2k + ll/4adu ) /(l - kr 2 ) 

- (ad + 2d 2 + 2k + ll/4adu ) r 2 

- (ad + 2d 2 + 2k + ll/4adw ) r 2 sin 2 6> 



(3.15) 



The geodesic deviation equation in the case of a perfect fluid along the neigh- 
boring world lines can be generalized within the Finslcrian framework (£ M is the 
deviation vector) [231 US] 



Within a Finslerian space-time framework the concept of constant curvature 
K is formulated bv[2^,[3%] 



4 Einstein's field equations for an anisotropic 
universe 

4.1 The energy- momentum tensor and the Friedman-like 
equation 

The energy- Momentum tensor of a Finslerian perfect fluid for a comoving observer (2 1( 
[35l [36] is defined to be 



where P = P{x), \x = ijl(x) is the pressure and the energy density of the cosmic 
fluid respectively. The vector y a = is the 4-velocity of the fluid since 
y a = (1,0,0,0) with respect to comoving coordinates. Thus T^ v becomes^,, = 
diag(p, -Pf tj ) in matrix form)[9l E3 [39l f4Tj 



where the operator 





(3.17) 



T^(x,y{x)) = + P)y,j,(x)y v (x) - Pf^(x, y{x)) 



(4.1) 




(4.2) 



T 



f°°H - 3P 



The substitution of (|4.1|) to the field equations 




(4.3) 



implies the following equations at the weak field limit 



G 



'^ + - A -uo = - ^(M + 3P) (4.4) 
a 4 a J 



a a k lid 
- + 2— + 2 — + — -u = 47tG/i-P 4.5 
a a z a z 4 a 



after subtracting (|4.4[) from (|4. 5[) we obtain the Friedman-like equation 



where we define Zt as 



a. \ 2 a 8irG k 

+ -z t = — f i-— (4.6) 
a 3 a z 



z t = ii . (4.7) 



The quantity z t is a constant since we consider a linear approach for the <fi{x) by 
(|2.5p . The previous equation is similar to the one derived from the Robertson- 
Walker metric in the Riemannian framework, apart from the extra term ^z t .We 
associate this extra term to the present Universe's anisotropy.In case we study 
Finslerian models with a cosmological constant ,44| the field equations (|4.4p . (|4.5p 
can be given in the following form 

- + —*> = - ^(/i + 3P) + A (4.8) 
a 4 a a J 

a a 2 k lid 
- + 2— + 2— + — - 

a a z a z 4 a 

and we end up with the equation of motion 



2^ + — -uo= 4irG((j,- P)+A (4.9) 



d\ 2 d 8ttG k A 



, . z t = —^-- + - (4.10) 
a J a 6 a z 6 



4.2 The parameter z% and the weak linearized anisotropy 

The physical quantity z t describes the variation of anisotropy which evolves lin- 
early due to (|2.5[) ;it depends upon the scalar (f>(x) which is the only quantity 
of the Lagrangian that gives us insight about the evolution of anisotropy. The 
parameter zt is measured by the Hubble's units as (|4. 6|) implies. It is signifi- 
cant that zt depends on the geometrical properties of the Finslerian space-time 
manifold. Indeed, the component Cqoo can be directly calculated from (|3.8[) as 



Good = y (4.11) 

and after differentiating with respect to proper time we lead to the direct de- 
pendence of z t on the Cartan torsion component Cqoo 

zt = 2C oo,o (4.12) 

hence the variation of anisotropy is closely related to the variation of the Cartan 
torsion tensor as an intrinsic object of the Finslerian space-time. 
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5 The cosmological anisotropic parameters 



We list the main anisotropic parameters constructed within the Finslerian framework |21[ 
[22l[28l[35] 



5.1 The anisotropic scale factor a(v(s)) and the Hubble 
parameter H 

The anisotropic scale factor is defined along each world line. S(s) is the length 
scale introduced in pj [5T] defined as 

S(s) = ~a(v(s)) (5.1) 

where v(s) is the tangent vector field along the world lines. The anisotropic 
Hubble parameter H is given by the relation(5 = ^■^(s)) 

B = § = i8 (5.2) 

The term is the expansion in the Finslerian space-time expressed as 

© = l& - C^p (5.3) 
The anisotropic Hubble parameter can be computed from ()4.6p 

= H 2 + Hzt (5.4) 

thus Hubble's units have to be attributed to Zt- Since H 2 > we should fix 
H 2 ^> \Hzt \ Therefore the lower limit of Zt is —H, 

z t > -H. (5.5) 



This result is in agreement with the fixing of a similar parameter in fill 1121 [13] 
for self- accelerated brane- world cosmology. Since the parameter z t is related to 
the variation of anisotropy we expect it to have negative sign(se// accelerating 
universe) as it may control a transition of the universe from a state of anisotropy 
to a smoother isotropic phase [IB]. 

5.2 The density and the deceleration parameter 

The density parameter can be defined with respect to the parameter H 

Cl u = — —a — - — (5.6) 

^ -Mi 2 Mcrzt ^ ' 



where 



U ~ 8ttG' 
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The deceleration parameter is defined in terms of the anisotropic scale factor 
S(s) 

q = -SS/S 2 . (5.8) 

The Friedman-like equation can also be rewritten in the form 

% - 1 = k/(H 2 a 2 ) (5.9) 

If /i < ficrit then f2 M < 1 or k < (open universe). If fi > fl cr it then f2 M > 1 
or k > (closed universe). The latter case /j, = /2 cr it corresponds to S1 M = 1 
or fc = f/!a£ universe) . Thus the values of influence the type of spatial 
curvature. In case we need to express H in terms of the redshift, the present 
value of the Hubble parameter H and the Qm ^a we insert [36l [38l [39] 



H(z) = H E(z) (5.10) 
into (|5.4[) where the quantity E(z) = H(z)/H is given by 



e(z) = y/nT t + V^z t + fiAfo (i + zf + n K „(i + z) 2 + n Ao (5.ii) 

and the parameter f2 Zt is defined by 



Therefore we are dealing with an expression for H (z) which depend up on the 
anisotropic parameter zt, the redshift z and the f2's. 

The Friedman equations in terms of Q 's 

The Friedman equations can be expressed in turns of the density parameters. The 
STs give back some useful insight about the unknown parameter z t and en- 
able us to test if our cosmological model fits to the current data(e.g.WMAP 
data). Indeed, the equations of motion (|4.6[) . (|4.10[) are reduced to the form 

1 = Q Mo + n Ko + fi Ao - z t /H . (5.13) 

In case we are interested in manipulating the Friedman equation for a specific 
value of the redshift z we can make use of (|5 . 1 1[) .It is difficult to distinguish 
the last term of (|5 . 13[) from GR dark energy, since both terms accelerate the 
universe. If the self-accelerating term z t H dictates over the expansion, the 
cosmological constant in (|5.13|) . (|5.11|) has to vanish [T2] . 

5.3 The continuity equation 

The energy density can be calculated by direct integration of the continuity 
equation = 0. A proper manipulation of (|4.1[) and (|3.9p leads to the equa- 
tion 

T° + (a/a + uo/2) (3T ° - 7?) = (5.14) 
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which can be simplified to the form(c = 1) 

/i + 3a/a(^ + P) +m /2(3P + 2^) = (5.15) 
since we have applied the approximations 

/ 00 = 2/(2 + u + 2 U 2)«l 

/°°= -2uo(l+4uo)/(2 + uo + 2u§) 2 

r b = / 00 A + /, oV«A-^o/2 [ ' 

3T ° — T\ w + 3P. 

A perfect fluid relevant to cosmology obeys the equation of state 

P = w[i (5.17) 

where w — for a matter dominated universe and «; = 1/3 for a radiation 
dominated universe. The substitution of (]5 . 1 7[) to (|5. 15[) leads to 

(i/n = -3(1 + w)a/a + ii /2{2 + 3w) (5.18) 

and the integration of the differential equation (|5. 18|) implies 

H oc _3(1+u,) exp {~u {3w + 2) /2) (5.19) 

therefore 



a 3 exp(— i*o ) matter dominated universe 
1 a -4 exp(— 3uo/2) radiation dominated universe 



(5.20) 



This asymptotic behavior indicates that the weak anisotropy affects homogene- 
ity. 

5.4 The look back time at the presence of weak anisotropy 

Given the definition of the anisotropic Hubble parameter (|5.4[) we can generalize 
the concept of the lookback time to — i* in 39J 

r z * dz' 

to-U= / —~ (5.21) 

Jo (l + z')H(z') 

where t$ is the age of the universe today and t* is the age of the universe when 
the redshift was z = z*. Therefore the general expression for the world time 
evolution should be [38 

f°° dz' 

t( z )= / = (5.22) 

since z — » oo at the start of the universe. Taking into account the weak anisotropy 
scenario and the redshift expression of the Hubble parameter (|5. 10[) into account 
we deduce the small zt expansion for the worldtime today to(z = 0) 

oo 

to = to + TkZ t ( 5 - 23 ) 

fc=l 
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where 



5.5 The Finslerian Raychaudhuri equation 

The generalization of Raychaudhuri's equation has been given by the following 
formula [21] 

<^> = -i<9 2 - a^a^ + Q^u>^ - AirGifi + 3P) + [y^y*)-^ (5.25) 

where a 2 — a^ v a^ v , Co 2 — uj^uj^ 1 ' are the Finslerian shear and vorticity respec- 
tively, which are defined in 33, 35Jand refer to a perfect fluid. The Eq. (|5.25p is a 
direct application of the Finslerian Lie derivative for dust-like matter developed 
in [24]. It includes the anisotropic gravitational influence of the matter along 
the worldlines which is expressed by the tidal force of the field 

= 4ttG(m + 3P). (5.26) 

Using the expressions (|5.8|) , (|5-3[) we produce the linearized Raychaudhuri equa- 
tion( q = — §-&~ 2 where H depends on a, a in virtue of (|5.4p ) 

\q 2 6 = 4ttG(m + 3P)^ + /(a, d, a, z t ) (5.27) 
3 H 

where the function /(a, d, a, z t ) is defined as 

f(a,a,a,z t )= -3H 2 - f § (tf 2 + ff 2 )+ 

+ i [27rG(/i + 3P) + 3z 4 #]z 4 . 

Thus we can expand for small values of z t and conclude that 

\e 2 q = 4*GQ* + BP) + 3% - ^a(t) - a(t) 2 + 
3 a 8 a (i) 

The sign of the right hand side of (|5.27|) determines the state of expansion. If the 
inequality /(a, d, a, z t ) < is valid then the term /(a, d, d, Zt) contributes to the 
acceleration of the universe(the field assist inflation), whereas if /(a, d, d, z t ) > 
it will slow the expansion down(the inflation domination must be longer to 
accelerate the universe). This specific effect is due to the kinematical reaction 
of the geometry of the spatial hypersurfaces, rather than an attempt to sup- 
press inflation(it can be considered as an essential ingredient of the Finslerian 
ansatz) [HI [33J • 
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6 Einstein field equations with anisotropic term 



We assume z t to be a constant and study the differential equations (|4.4p . (|4.5p 
and (14. 6|) both for the cases of matter and radiation dominated universe. We 
notice that for z t = the field equations reduce to the usual ones coming from 
a Riemannian Robertson Walker metric [371 SD S3]- The whole calculation is 
done for a homogeneous universe of constant density [i. 



6.1 Solution for a matter dominated universe 

The solution for the scale factor in the case of a matter dominated universe 
is derived by the integration of the Friedman-like equation (|4.6p with initial 
condition 

o(0)=0 (6.1) 

where z t = const. We set up t = as the beginning of time without considering 
any quantum effects;there are different ways of handling the initial condition 
for example setting a(tpi) = or considering the scale factor after the Plank 
scale(e.g. see [30]. The Eq. (|4.6|) is simplified if we insert the parameter (confor- 
mal proper time) j37j [4TJ 

5= H (6-2) 
Jo a(ui) 

which measures the arc in rad traveled along by a photon on a sphere of radius 
a(t).We study the asymptotics of the scale factor 

a(t) = a(t(r})) = a(n) -> (6.3) 

where cdt — adv. We consider a homogeneous universe of constant density /i 
calculated as [JT] 



M M 
V ~ 2iT 2 a 3 

We assume that the universe would take over the same volume as in the closed 
case (k = 1) since (I5.20P is valid for the asymptotics of u. The parameter 
k determines the kind of spatial geometry k = 0, — 1, +1 (flat, open, closed 
universe respectively). Since - — a~ 2 a Eq. (|4.6p becomes 

a + -fa 2 = 4« ~ ^ + ^ — a 6.5) 
V 2 / 4 37T 

We study a Universe that accelerates very fast at its early stages thus we can 
accept a+ya 2 > 0. The velocity of expansion a takes on very large values. The 
Eq. (|6.5|) can be integrated directly for all the values of k. It is more convenient 
to substitute [37] 

a * = ~3^r ^ ' 
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and the separable Eq. (|6.5p leads to 



2 1 .L kx- 2a* 



together with the initial condition 

a(0) = 0. 



-dx 



I) calculation for k = 

We fix fc = at JO) and find 



6ff * 



Io = Jo" V^ 4 + &a*/z?xdx 
we expand for a small Zt and arrive at the solution 
V2 



t = 



3-,/a" 



: a 3/2 + l/(12a*)a 3 z t /c + 0(z 2 ) 



II) calculation for k = — 1 

A direct integration of (!6.7() leads to 



t = ~ [a 2 /2 + 2a*a + 4a 2 log |a - 2a*| + Ji] 



where 



_ r y/z 4 - 4/z 2 x 2 + 8a*/z 2 ~r 
Jo x 2a* 



after expanding for small Zt 



t= y/ a(a + 2a*) — a* log ^1 + a/a* + \J a(a + 2a*)/c 
- (a 2 /4 + a*a + 2a 2 log |a - 2a* |) z t + O(zf). 



Ill) calculation for k = +1 

The calculation is the same as the previous case 



t = y [a 2 /2 - 2a*a + 4a 2 log \a + 2a*| + J_i] 



where 



_ r y/x 4 + 4/z 2 x 2 + 8a„/zfx 
1 Jo x + 2a* 
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and after expanding for small z t we obtain the solution 



— y/a(2a* — a) + a* arccos (1 — a/a*) + 
(a 2 /4 - a*a + 2a 2 log|a + 2a*!) z t + 



(6.16) 



The leading term of the solutions for all k represents the solution given by 
the field's equations of the Robertson- Walker metric [57]. A small a expansion 
gives the asymptotic behavior t ~ a 3 / 2 or equivalently 



,2/3 



(6.17) 



6.2 Solution for a radiation dominated universe 

The solution for a radiation dominated universe can be deduced by inserting 
the equation of state P = \^ ra d into (|4.4p and after adding this to f|4. 5|) we end 
up with the equation 

d T 

— {ah) = —k — -aazt (6.18) 
which we integrate and find 

7 

da+ -z t a 2 = -kt + Ci (6.19) 
4 

If we substitute z t = to (|6.19p we get back the usual solution for a radiation 
dominated universe a oc y/t for all values of /c.If we take into account the initial 
condition a(0) = we arrive at the solution 

a(t) = |(C 2t +k)(l — exp(-^z t t)^j - \kz t t^ ' . (6.20) 
The expansion of the solution for small z t is 



6.3 Solution for the de-Sitter model 



(6.21) 



The de-Sitter model for an empty anisotropic universe constructed in|21j leads 
to the equation of motion(/z = in (|4.10p ) 



a\ 2 d k A 



- ' T7 (6-22) 

Since H 2 in (I6.22p can be written as H 2 = (^) 2 + > 0, the cosmological 
constant A is restricted by the inequality A > 3K where K = k/a 2 is the 
curvature of the space. 
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I) calculation for k = 

The case of zero curvature can be integrated to give 



1 /2 

i(t) = const x exp (-z t /2t) exp (A/3 + z\jA) t 



(6.23) 



The solution converges to the one without anisotropy if we let z t — > 0|44j. 
II) calculation for the special case A = 3K 

The field equations for a space of constant spatial curvature(maximal symmetry) 
and an empty universe(T Alt/ = 0) imply the condition[2Tl [39] 

A = 3K (6.24) 

hence A can be substituted to the equation of motion (|6 . 22[) 



. \ 2 

- + -z t = 6.25 
a J a 

and the solution for the scale factor is 

a(t) = const x exp (—z t t) . (6.26) 

In both cases the constant of integration can be absorbed into a(t) if we choose 
the right scale(e.g.o(0) = 1). 

6.4 A model for inflation with anisotropy 

The idea of inflation can be incorporated into the model of Friedman-like equa- 
tions with weak anisotropy if we introduce the vacuum energy density of a scalar 
field Vq to the energy density fi. Indeed, we consider Gu — Vo/ TO p;> i-e. £t = Vq 
thus (|MD implies [31 SD1I11] 

• 2 , • 8"7T Vq 2 

a + aaz t = 5— a — k. (6-27) 

3 m z pl 

Since we work at an inflationary phase the size of the scale factor is such that 
the term ^--^-a 2 dominates over k hence we can neglect k and (|6.27[) can be 
rewritten as 



H 2 = Y ^- ( 6 - 28 ) 

o rn pl 

Taking the positive square root of (|6 . 28[) we find the scale factor 

a(t) = const x exp (—z t /2t) exp (87r/3 • Vo/mpi + zf /4) 1/2 ij (6.29) 

recovering the de-Sitter solution (|6.23|) and the expected exponential rate of 
expansion for the early inflationary phase of the universe. 
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7 Estimation of the Cosmic Microwave Background 
Radiation (CMB) 

The estimation of CMB can be achieved with the aid of Stefan-Boltzmann's 
law [SUSS] 

MradC 2 = ctsb9*T(z) 4 (7.1) 

where 

7T 2 k 4 B . . 

aSB = 30^' (7 ' 2) 
The temperature T(z) is the radiation temperature for a given redshift z and 

bosons fermions 

is defined as the sum of the boson and fermion spin states(e.g. for photons g = 2, 
for neutrinos 3 = 1 and for massive particles g = 2s + 1 ).The calculation could 
be attainable if we consider some data for the Hubble parameter H = a/a and 
the anisotropic constant z t .The radiation dominated solution does not depend 
on the nature of the spatial geometry of the universe thus we fix k — in (|4.6|l 
and obtain fj, ra d 



Vrad - $^(H 2 + Hz t) 
therefore l[7T|l .([7^] l yield [37] 



a' fj2 (7-4) 



1 / 3 



The calculation can be directly derived if we manipulate (|7.5p into the form 

k 2 B T{zf = y L H M (7.6) 



where mpj stands for the Planck mass 

Thus we can calculate the temperature T(z) for a given value of the redshift z 
with the aid of the formula 

T(z)=T H (z)[l + z t /H(z)} 1/4 (7.8) 
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where Th(z) is the value of the temperature for a given value of the Hub- 
ble parameter without the assumption of weak anisotropy and H(z) is calcu- 
lated from (l5.10p .Due to the adiabatic expansion of the universe we can trans- 
form the value T(z) at the redshift of CMB to its present value To using the 
formula(a C MS = 1/(1 + z R ), z R w 1090) 

T = clcmbT{z). (7.9) 

8 Discussion 

The study of a FRW-model with a weak vector field incorporated in the met- 
ric structure of space-time provides us the extended Friedman- like Eq. (|4.6[) .Thc 
contribution of the variation of anisotropy is expressed by the additional pa- 
rameter z t produced by the Finslerian character of the geometry of space-time. 
Especially as it is evident from Eq. (|4.12[) z t has a direct dependence up on the 
Cartan torsion component Cooo- We remark that our present model correspond 
to the ones studied in [T3] H3] for a hat universe due to the correspondence of Zt 
to where r c is the extra parameter defined there. The extra parameter z t 
appears to compete against the contribution of the cosmological constant due 
to (|5TT5|) . 

We perform the model- independent and insensitive to perturbation S'-test[4], 
where S is the CMB shift parameter p~2l [13] .We are pretty confident that our 
model reproduces the WMAP data for a flat universe since for S = 1.70 ± 
0.03 [16] we end up with \SIjc\ <C l.The same result seems to be valid if we 
apply the baryon acoustic oscillation peak test for A = 0.469 ± 0.017(5], [3]. The 
procedures of the tests and the formulas for A and S can be found in [12]. A 
part of a future work is the investigation of cosmological perturbations since 
the Finslerian approach generates deviations from homogeneity and isotropy. 
A more fundamental task is the comparison of our model to the data of CMB 
anisotropics and the matter power spectrum, strongly connected to the analysis 
of the density perturbations [14] . 

The initial highly compressed, thermal radiation dominated state of anisotro- 
py is considered to be adiabatically transformed to a cooler matter dominated 
isotropic phase [41] .In such a case where the anisotropy energy is converted to 
thermal energy and large amount of entropy [46], a phase transition in the ge- 
ometric structure can be regulated by the second term of Eq. (|2.ip . since we 
expect 4>{x) monotonically decrease as the universe expands, in order to obtain 
the standard FRW model. This ensures negative values for the parameter zt 
and provides us a self accelerating cosmological model. 

The whole picture of anisotropy directed by a primordial vector field can be 
locally incorporated to the anisotropic metric structure of a Finslerian space- 
time. The osculation of the Finsler space leads to the construction of a disformal 
Riemann structure and can be interpreted as a model of modified gravity. 
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A Appendix 

In the following we present some basic elements of Finsler geometry [l8l [T9l [26l 
[2Z1I33]- In 1854 B.Riemann, before arriving at Riemannian metric was concerned 
with the concept of a more generalized metric 

ds 2 = F(x 1 ,x 2 ,...,x n ,dx 1 ,...,dx n ) (A.l) 

where n is the dimension of the space. A Finsler structure is provided by a 
n-dimensional C°° manifold M n , a C°° function F = F(x,y) defined on the 
tangent bundle T~M = TM/{0}, F : T~M -> R that satisfies the conditions 

(Fl) f(i,y)>0V^0 

(F2) T(x,py) = pJ r (x,y) for anyj? > ^ ' ' 

where y denotes the directions or velocities on the considered manifold with the 
previous coordinates. The metric tensor (Hessian) 



1 d 2 T 2 

2 dyWyJ 



fijfav) = o a ..ia.A x >y) ( A - 3 ) 



is of rank[(fij)ij] = n and homogeneous of zero degree with respect to y due 
to the Euler's theorem. The length s of a curve C : x l (t),a < t < b on the 
manifold is 

b 

T{x(t),y{t))dt. (A.4) 

The integral of the length is independent of the parameter if and only if the 
condition (F2) is valid. The condition of homogeneity enables us to define the 
line element 

ds = T(x, dx) (A. 5) 

and the variation of the arclength 8 J ds = implies the Euler-Lagrange equa- 
tions — = which represent the geodesies of the Finsler 
space. The equation of geodesies then becomes analogous to the ones of the 
Riemann space 

d 2 x i 



ds 2 



lWv k = (A.6) 
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where the Christoffcl symbols are defined by the usual formula 

ijk( x >v) = \r r { x ,y) (frjA x >y) + frk,j(x,y) - fjkA x ,y)) ■ ( A - 7 ) 

The notion of torsion tensor is crucial within the Finsler Geometry's frame- 
work. A Finsler space is a Riemann space if and only if C'ijk — where Cijk is 
the torsion tensor defined by E.Cartan as 



Cijk (A " 



Therefore a Finsler space can be treated as a natural generalization of a Riemann 
space. 

B Appendix 

For a Finslerian vector field X a (x,y(x)) the S— covariant derivative has the 
formpH [35] 

X* p {x,y{x)) = X%{x) + T™{x,y{x))X>>{x) (B.l) 
and the Cartan's covariant derivative is given by 



Xfo{x,y{x)) = X%x) - ^(x,y(x))G p (x,y(x))+ 
+ T^(x,y(x))XP(x). 



(B.2) 



The are the Cartan's connection components defined as 

nU x > y) = - C%G% - C^G P X + C Xw GPg VK ) {x, y) (B.3) 
and the G»,G» are 

Gt = (B.4) 

2G" = 7 ^y p y CT - (B.5) 
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Zt/H 


Q. K = 1 - M + z t /H 


0.110 


-0.9029492854 


-0.0129492854 


0.151 


-0.8443621881 


0.0046378119 


0.232 


-0.7539951109 


0.0140048891 


0.244 


-0.7710952485 


-0.0150952485 


0.251 


-0.7608626761 


-0.0118626761 


0.252 


-0.7327905974 


0.0152094026 


0.253 


-0.7610339645 


-0.0140339645 


0.254 


-0.7556133186 


-0.0096133186 


0.255 


-0.7449999998 


0.00000000002 


0.259 


-0.7313813921 


0.096186079 


0.260 


-0.7400000040 


-0.000000004 


0.261 


-0.7236653148 


0.0153346852 


0.262 


-0.7378258680 


0.0001741320 


0.263 


-0.7464096424 


-0.0094096424 


0.264 


-0.7249493659 


0.0110506341 


0.2651 


-0.7441489179 


-0.0092489179 


0.270 


-0.7425851580 


-0.012585158 


0.271 


-0.7389289128 


-0.099289128 


0.272 


-0.7279999954 


0.0000000046 


0.273 


-0.7337402645 


-0.0067402645 


0.274 


-0.7135079061 


0.0040405097 


0.275 


-0.7151439555 


0.098560445 


0.2761 


-0.7238999999 


0.0000000001 


0.277 


-0.7363820256 


-0.0133820256 


0.278 


-0.7166084631 


0.0053915369 


0.279 


-0.7339876359 


-0.0129876359 


0.281 


-0.7038915738 


0.0151084262 


0.282 


-0.7096732356 


0.0083267644 


0.2831 


-0.7079203819 


0.089796181 


0.284 


-0.7285702833 


-0.0125702833 


0.285 


-0.7232538088 


-0.083538088 


0.2861 


-0.7229838341 


-0.090838341 


0.289 


-0.710999993 


0.000000007 
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0.290 


-0.7032528013 


0.0067471987 


0.291 


-0.7039021795 


0.0050978205 


0.295 


-0.7002433383 


0.0047566617 


0.2981 


-0.6924644204 


0.0094355796 


0.325 


-0.6674671885 


0.0075328115 


0.341 


-0.6560099554 


0.0029900446 


0.345 


-0.6501638547 


0.0048361453 


0.351 


-0.6617634475 


-0.0127634475 


0.36 


-0.6545617924 


-0.0145617924 


0.372 


-0.6385195052 


-0.0105195052 


0.38 


-0.6199999969 


0.0000000031 


0.392 


-0.6008829186 


0.0071170814 


0.411 


-0.5985044824 


-0.0095044824 


0.421 


-0.5789999998 


0.0000000002 


0.451 


-0.5563735170 


-0.0073735170 


0.512 


-0.4989395368 


-0.0109395368 



Table 1: The A— test estimation for the parameter Zt-We insert the values of the 

r H 3 d 2 Iz-,) "I 1 / 3 

Q.M to the baryon acoustic oscillation peak A = v"M H z °^(i^_ zi )2 an£ l calculate 
back the corresponding values of zt, where A = 0.469 ± 0.017 3 , Z\ = 0.35 the typical 
luminous red galaxies redshift and cLl(z) the luminosity distance defined in various 
Relativity textbooks(e.g. see [35]). The parameter Qk close to zero as expected from 
the present observational predictions for a flat universe of WMAP [16| . 
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Figure 1: The zt parameter versus flM;the A-test values of —Zt/Ho are depicted 
by the "plus sign" line and compared to the "dashed line" values coming from the S- 
test.Both tests almost reproduce the theoretical values of z t for a flat universe(f2A' = 0) 
calculated by the Friedman equation —z t /Ho = 1 — SI m- ) 
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